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A method of analyzing, in detail, the "performance of symmetrical Cas- 
segrain antennas has been developed that uses a digital computer efficiently. 
For a specified antenna geometry and feed excitation, the program will 
compute and graphically display the amplitude and phase illumination 
of the subreflector, main reflector, and far-field pattern. These results may 
be used to optimize antenna performance by changing parameters and ob- 
serving the effect. 

Analysis of a Cassegrain antenna with a near-field conical horn feed is 
discussed as an application of the method. Because the radiation character- 
istics of the horn are determined by the horn flare angle rather than the 
horn aperture, broadband performance is obtained. It was indeed found that 
a 50 per cent bandwidth is achieved with a dual mode TE lx — TM 1X mode 
feed, provided the proper phase relationship between the modes can be 
maintained over the band. For dual mode excitation an aperture efficiency 
of 70% and a noise temperature due to the power loss at the sub and main 
reflectors of less than 6.5°K was obtained. For a single mode feed (TE^), 
there was a degradation in the E-plane side lobe levels and a corresponding 
10°K increase in noise temperature. Excitation in the TM 01 mode was also 
examined for angle-error sensing purposes. Also, the antenna can be used 
with reasonable efficiency well below the design frequency in which case it 
functions as a far-field fed Cassegrain antenna. 

I. INTRODUCTION 

The essential radiation characteristics of multiple reflector anten- 
nas can be predicted very accurately with existing analytical and 
computational methods. Previous work on the open Cassegrain antenna 
showed that good agreement can be achieved between calculated and 
experimental results. 1 Deviations occurred mainly in the sidelobe 
regions of the radiation patterns, and these had only a small effect on 
overall antenna performance. 
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We are concerned here with a simpler problem, but one of perhaps 
more general interest: the analysis of symmetrical Cassegrain anten- 
nas. We present a computational method in which the amplitude and 
phase illuminations of the subrefiector and main-reflector, as well as 
the far-field radiation pattern, are determined in detail, given the 
geometry of the antenna, the dimensions of the feed horn, and the 
excitation modes of the feed. The analysis includes near-field excita- 
tion — an important configuration for broadband operation. Included 
in the program is a graphic routine which plots all radiation patterns, 
the intermediate illuminations and the final far-field results. Because 
of the this feature, and the fact that only seven parameters arc 
required to define the geometry of the antenna, the program is par- 
ticularly useful for optimizing antenna performance. 

The symmetry of the antenna results in improved computational 
efficiency. For the open Cassegrain antenna, double-integration was 
required to compute radiation patterns. An approximation recently 
was obtained 2 which, when applied to symmetrical Cassegrain anten- 
nas, eliminates one integration with only a small reduction in ac- 
curacy. 2 This makes it possible to compute the radiation char- 
acteristics of large Cassegrain antennas a few hundred wavelengths 
diameter in minutes. 

II. NEAR-FIELD SYMMETRIC CASSEGRAIN ANTENNA 

The antenna under consideration was intended to be used as the 
ground station of a satellite communications system. Wide-bandwidth 
(25 per cent) and low-noise requirements motivated the choice of a 
near-field conical-horn symmetric Cassegrain configuration. The 
near-field feed produces relatively low spillover at the subrefiector, 
resulting in a lower noise temperature. 3 Also, because radiation of the 
feed is virtually confined to the geometrical illumination region of the 
horn, 4 there is a larger potential bandwidth available. 

An additional requirement had to be explored: operation at about 
!/4 nominal frequency, for target-acquisition, by using both TE n and 
TMoi mode excitation. This leads to the choice of a near-field design 
at the higher frequency because it would tend to function as a con- 
ventional far-field design at the lower frequency. At the lower fre- 
quency there would, however, be a shift of the phase center towards 
the horn aperture, resulting in a phase error in the subrefiector il- 
lumination and a consequent reduction in efficiency, but perhaps it 
would be adequate for the intended function. Finally, dual-mode 
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illumination using the TE n and TM n modes was of interest because 
of the nearly circular symmetric radiation patterns that can be 
obtained. 5 

III. CASSEGRAIN ANTENNA GEOMETRY 

Figure 1 shows the geometry of a Cassegrain antenna. It consists 
of a conical feed horn, a hyperboloid subreflector and a paraboloid 
main reflector. One focal point of the hyperboloid coincides with 
the focal point of the paraboloid and the other focal point with the 
phase center of the horn. The main reflector illumination angle is 
equal to the geometrical subreflector illumination angle, 8 m . The feed 
is located in the geometrical shadow region of the subreflector. 

The initial design of a Cassegrain antenna is usually based on 
geometrical optics, which imposes certain restrictions on the antenna 
geometry. The constraints are that the feed horn be located in the 
shadow region of the subreflector and that the subreflector intercept 
most of the power radiated by the horn. 

To relate the radiation properties of the horn to the antenna 
geometry it is convenient to define a parameter K by: 

K = ^smS (1) 

where 

d = horn aperture diameter 

X = wavelength 

8 = the angle subtended by the subreflector with respect to the 
center of the horn aperture (Fig. 1). 

For conventional Cassegrain antennas representative values of K 
are from 1.2 to 1.6. For these values of K the major portion of the 
main lobe of a narrow angle horn excited by TEn and TMn modes, 
is intercepted by the subreflector. The lower value of K is preferable 
for TEn mode excitation because the major lobe of the horn radia- 
tion pattern is narrower in the E plane than in the H plane. Beyond 
the major lobe region the phase variations are too large for efficient 
subreflector illumination. 

For near field Cassegrains the values for K are much larger, such 
that the radiation characteristics of the horn are primarily determined 
by the horn flare angle. 

For the feed horn to be located in the geometrical shadow region of 
the subreflector it is necessary that angle 6 b be not less than angle 6 bh , 
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Fig. 1 — Geometry of cassegrain antenna. 

that is, 

tan 6 b = B tan 6 bh (2) 

where B is a constant with B > 1. Figure 1 shows the angles 6 b and 
dbh- The horn blocking angle 6 b u can be expressed in terms of the 
geometrical parameters and (1) by 



tan d hh = 



K\ sin 6 r 



(3) 



7bh D sin (5 + O 
With all other parameters specified, d bh is minimum for 

8 + m - |- (4) 

Similarly, with 6A also specified, if has a maximum value when (4) holds. 
Using (2) and (3) and expressing d b in terms of the geometrical 
parameters, the following equation is obtained for the subreflector 
diameter, D. 



D = 



2BKf\ sin g. 



'sin(0 m + 8) +^fram <>,. 



(5) 



16/ 
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where 

/ = focal length of the paraboloid main reflector. 

Equation (5) agrees with the previously given condition for no block- 
ing 6 D zz (2Kf\) Vi . In practical antenna designs BK\/lQf is small 
compared with sin (8 -1- d,„), hence (5) may be rewritten in terms of 
the main reflector diameter, D , as : 



D > cos yj_M_ (6) 

\2/\ s in(5+0 

Equation (6) shows that antennas with large main reflectors also 
require larger subreflectors, but that the ratio (D/D ) z which is a 
measure of the amount of power blocked by the subreflector is in- 
versely proportional to the main aperture diameter D . Hence the 
condition (6) is of importance primarily in the design of relatively 
small Cassegrain antennas. 

Equation (6) also shows, as expected, that a conventional Casse- 
grain antenna requires a smaller subreflector than a near-field Cas- 
segrain antenna, since K is smaller for the former. However, this dis- 
advantage of the near-field Cassegrain antenna is offset by other 
advantageous properties. 

Another parameter which influences the antenna design is the 
total Fresnel number of the horn at the subreflector distance, defined 

by 

'.-Sft+fr 

For a conventional Cassegrain, F t can be selected, to a certain ex- 
tent, independently of the antenna geometry, because the radiation 
properties of the horn are not directly related to the horn length, I. 
For such an antenna with combined TE U and TM n mode excitation, 
a total Fresnel number in the 0.5-0.65 range would provide a nearly 
uniform subreflector illumination over a wide frequency range (about 
30 per cent) with relatively small phase deviations. For TEn mode 
excitation, a lower Fresnel number is necessary, because the phase 
of the S-plane horn radiation pattern is more frequency sensitive for 
larger Fresnel numbers. 

For a near field Cassegrain antenna the total Fresnel number is 
almost directly related to the antenna geometry. Specifically, for an 
antenna with the horn located in the shadow region of the subreflector 
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and with a subreflector illumination angle equal to the horn flare 
angle, Ft is given by 



with 



/, + I. 



tan 6 b = ?-r-- (9) 



®" 



1 w 

Since Z is much smaller than k the total Fresnel number, F t , is 
mainly determined by the subreflector diameter D. 

For a near field Cassegrain it is necessary to have both K and F, 
large. Equations (3) and (8) show that these quantities are propor- 
tional to D 2 . 

IV. ANTENNA DIMENSION 

For the antenna under consideration the main reflector dimensions 
were specified. Its diameter, D is 224a (A = wavelength at the design 
frequency), its focal length, /, is 72.8a and the corresponding geo- 
metrical illumination angle, 6 m , is 75°. 

The initial choice of the other antenna dimensions was based on 
the following consideration. As shown above, K has a maximum for 
8 + 6 m = ir/2. Using this condition the subreflector diameter, D, has 
been chosen such that the optimum value of K is unity at the lowest 
frequency (\ L = 4.5A). At the design frequency, K is about 3 times 
larger than is required for a conventional Cassegrain antenna feed. 
For this value of K, D is 25a. With these parameters a horn with a 
maximum diameter of 17.6A can be located in the shadow region of 
the subreflectors. The corresponding horn length, I, is 100A. However, 
a feed horn with these dimensions would introduce, at the lowest 
frequency, appreciable phase variations at the subreflector owing to 
the shift of the phase center of the horn radiation pattern. For this 
reason these horn dimensions were not used in the computations. 

The horn dimensions used were d = 14A and I = 42.5a. With these 
horn dimensions K is only slightly less than the optimum value. The 
location of the phase center, which is 5A in the front of the horn 
vertex, and the subreflector illumination angle, which is less than 
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the horn flare angle, were chosen on the basis of the computed horn 
radiation patterns for combined TE al and TM n mode excitations. 
Table I summarizes the antenna dimensions. 

V. PROGRAM FOR COMPUTING ANTENNA CHARACTERISTICS 

Programs have been developed which compute the antenna radia- 
tion characteristics and plot the computed radiation patterns. The 
computational methods are similar to those used in the computation 
of characteristics of the open Cassegrain antenna. 1 However, only- 
single integrations are used; one integration was eliminated by using 
the Fresnel region approximation for wide angles and large Fresnel 
numbers. 2 Appendix A gives the equations used. Appendix B discusses 
operational aspects of the programs and gives flow diagrams. 

The antenna characteristics for combined TE n -TM u and TM , 
mode excitations are computed in one operation. The computer pro- 
gram consists of three parts which compute (i) the horn radiation 
patterns, (u) the subreflector radiation patterns, and (mi) the far 
field radiation patterns. 

The horn radiation patterns are computed at a constant radius, s, 
corresponding to the subreflector distance. From these computations 
the power loss at the subreflector is obtained by integration. The horn 
radiation patterns are also computed at the subreflector surface to 
obtain the subreflector illumination. 

The subreflector radiation patterns are computed at a constant 
radius, /, and at the main reflector surface. From these computations 
the power loss at the main reflector and the main reflector illumina- 
tion is obtained. 

From the computed main reflector illumination the aperture gain, 
aperture efficiency and finally the far field radiation patterns are 
obtained. 

The antenna gain and antenna efficiency are determined from the 

Table I — Antenna Dimensions 

Main reflector diameter, D 224X 

Focal length, / 72.8X 

Main reflector illumination angle, 6 m 75° 

Subreflector diameter, D 25X 

Subreflector illumination angle, 0a 9.5° 

Horn length, I 42 . 5X 

Horn flare angle, a 9.5° 

Phase center location, p 5. OX 
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computed aperture gain and efficiency respectively by including the 
loss at the subreflector and main reflector. 

In the computations the phase variations at the sub- and main 
reflectors are included. Also included is the effect of the main reflec- 
tor aperture blocking by the subreflector but not the effect of the 
protruding horn. 

An estimate of the antenna noise temperature is obtained by as- 
suming, somewhat arbitrarily, that near the horizon half the power 
lost at the sub- and main reflectors contributes to noise. At zenith 
it is assumed that the power lost at the main reflector contributes 
to noise. A ground temperature of 300°K is used in the computations. 
The additional noise from possible scattering of the subreflector sup- 
port and the noise from the wide angle sidelobes of the far field radia- 
tion patterns are not included in the computations. 

VI. COMPUTED ANTENNA CHARACTERISTICS 

The antenna characteristics have been computed with the above 
computer program for the following feed horn excitations: (i) TEn 
and TMn at the design frequency, /„, 0.8 f and 1.3 /„, and (ii) TE M 
and TMoi modes at /„ and 0.22 f . The antenna characteristics for the 
different modes and frequencies are summarized in Table II. The 
tabulated power losses are normalized with respect to the total power 
radiated by the feed horn. 

6.1 TE ii and TM lt Mode Excitations 

The computations were performed by assuming that the two modes 
are in phase at the horn aperture. The TMn to TEu power ratio 
was assumed to be 0.17. This value was used because at the design 
frequency it minimizes the phase variations of the horn radiation 
patterns at the subreflector both in the E and H planes. 

Computations have been performed at the design frequency, f , 0.8 
/„, and 1.3 f . This corresponds to about a 50 per cent bandwidth. 
Except for the expected change in the antenna gain, the antenna 
radiation characteristics remain virtually the same across this fre- 
quency range. This indicates that if a frequency-insensitive conical 
feed-horn using TEn and TMn mode excitation could be developed, 
this antenna design is capable of efficient radiation over a 50 per cent 
bandwidth. 

Figures 2 through 6 show the antenna radiation patterns at the 
design frequency, /„. Included are: the amplitude and phase of the 
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Table II — Calculated Antenna Characteristics 



Frequency 


Design, /o 


O.8/0 


1.3/o 


/o 


0.222 /■> 


Mode 


TE„ 
and 
TM„ 


TMoi 


TEn and TMn 


TEn 


TE„ 


TMoi 


Per cent of: 

Power loss at sub- 
reflector 


2.7 


20.0 


3.4 


2.1 


9.3 


26.0 


50.8 


Power loss at main 
reflector 


0.7 


2.9 


0.9 


0.5 


1.0 


4.4 


5.6 


Power blocked by sub- 
reflector 


5.0 


— 


4.8 


5.5 


3.3 


2.9 


— 


Aperture efficiency 


72.7 


— 


73.4 


71.5 


75.4 


82.2 


— 


Antenna efficiency 


70.4 


— 


70.3 


69.8 


67.7 


57.2 


— 


Antenna gain, dB 


55.4 


49.4 
at max 


53.4 


57.6 


55.2 


41.4 


33.3 
at max 


Antenna noise 

temperature, 

°K 


Near 
horizon 


5.1 


34.35 


6.45 


3.9 


15.45 


45.6 


81.6 


At 
zenith 


2.1 


8.7 


2.7 


1.5 


3.0 


13.2 


16.8 


First 
sidelobe, dB 


H plane 


-22.7 


-12.9 
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-24.1 


-17.2 
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Fig. 2 — Subreflector illumination, TEn and TMu modes, freq. = /„. 
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Fig. 3 — Horn radiation pattern at distance S (TEn and TMu modes, freq. 
= /.) 

subreflector illumination, the amplitude and phase of the horn radia- 
tion pattern at the subreflector distance, the amplitude and phase 
of the main reflector illumination, the amplitude and phase of the 
subreflector radiation pattern at the focal distance, and the far field 
pattern. 

These figures show that the phase variations of the sub- and main 
reflector illuminations are relatively small at this frequency. This is 
because the location of the phase center and the ratio of the TM n 
to TEn modes has been chosen to minimize the phase variations at 
the subreflector. These figures also show that the far field radiation 
pattern is virtually the same in the E and H planes. This should 
result in a nearly circular symmetric far field radiation pattern. 

Figures 7 and 8, and Figures 9 and 10 show some of the antenna 
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Fig. 4 — Main reflector illumination (TE U and TMu modes, freq. = /<,). 
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Fig. 5 — Subreflector radiation pattern at focal distance (TEn and TMn 
modes, freq. = /„). 

radiation characteristics at 0.8 f and at 1.3 f , respectively. The phase 
variations of the sub- and main reflector illuminations, though small, 
are larger than at the design frequency. This primarily results from 
the shift in the phase center of the horn radiation pattern. It is the 
shift in the phase center which ultimately limits the upper frequency 
of operation for this type of antenna. 

6.2 TM i Mode Excitation 

The radiation characteristics for the TM m mode excitation were 
computed at the design frequency, f . Figures 11 through 13 show 
representative radiation patterns for this mode. Particularly pro- 
nounced are the amplitude oscillations of the main and subreflector 
illuminations. This seems to be characteristic for the TM 0] radiation 
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Fig. 6 — Far field radiation pattern (TEi, and TM„ modes, freq. = 
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Fig. 7 — Subreflector illumination (TEn and TMu modes, freq. = 0.80 X /<>). 

patterns from apertures and reflectors which are large compared with 
the wavelength and which are illuminated with nearly spherical wave 
fronts. However, no experimental evidence has been found to confirm 
these characteristics. 

The advantage of this antenna for TM i mode excitation com- 
pared with a conventional Cassegrain is less spillover at the subre- 
flector, hence, less antenna noise. However, the sidelobe levels of the 
far field radiation pattern are perhaps a few dB higher than could 
be obtained with a conventional Cassegrain. 

6.3 TE U Mode Excitation 

In view of the difficulties in realizing a conical feed horn with TEn 
and TMn mode excitation which would maintain the proper phase 
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Fig. 8— Far field radiation pattern (TEu and TMw modes, freq. = 0.80 x /•) 
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Fig. 9 — Subreflector illumination (TE„ and TM U modes, freq. = 1.30 X /■>)• 

relationship at the horn aperture over a wide frequency range, TEn 
mode excitation only has been investigated for the same antenna 
geometry. 

Figures 14 and 15 show some of the radiation patterns for this 
mode at the design frequency, / . The computations show that the 
phase variations of the sub- and main reflector illuminations are 
considerably larger, particularly in the E plane, compared with those 
obtained by using combined TE n and TMn mode excitations. Also, 
the sidelobe levels of the far field radiation pattern in the E plane are 
considerably higher than in the H plane. 

Table II shows that the computed antenna gain is 0.2 dB lower 
than the computed gain for TE U and TM n modes. However, the 
most significant difference is the increase in the antenna noise tem- 
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Fig. 10 — Far field radiation pattern (TE U and TM U modes, freq. = 1.30 X /..). 
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Fig. 11 — Horn radiation pattern at distance S (TMoi mode, freq. = /«). 

perature by 10°K near the horizon. This increase is primarily caused 
by the larger power loss at the subreflector because of the E plane 
horn radiation pattern characteristics. 

A reduction of the antenna noise temperature by a few degrees 
might be possible by increasing the subreflector illumination angle 
perhaps even beyond the geometrical illumination angle of the horn. 
Figure 14 shows that the phase variations in the E-plane radiation 
pattern are not very large in the vicinity of the presently-used sub- 
reflector illumination angle of 9.5 degrees. The computed horn power 
radiation patterns show that if in the present design the illumination 
angle were 10.5 degrees the antenna noise temperature would be re- 
duced by 4.8°K. The antenna gain for such a design would be reduced 
by only a small amount. 
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Fig. 12 — Main reflector illumination (TMoi mode, freq. = /«). 



ANALYZING ANTENNAS BY COMPUTER 



911 



2 -18 



-24 

-30 

-36 
-42 
■48 
-54 

-60 



A 
















/ 


1 
















^ 
















\ 
















1 


V 


\ 














\/ 


\/« 














/ 


V 


\f\ 




\ 








f 


V 


/ 


/V 


\ 












y 


f 





















0.4 0.8 1.2 1.6 2.0 2.4 2.8 
ANGLE OFF AXIS IN DEGREES 



3.2 



Fig. 13 — Far field radiation pattern (TMn, mode, freq. = /»). 

The bandwidth characteristics for this mode in the vicinity of the 
design frequency should be similar to those of the combined TE n 
and TM n modes. 

6.4 TE n and TM 01 Mode Excitation at 0.22 f 

The horn and far field radiation patterns for these modes are shown 
in Figs. 16 through 19. The right side of Table II summarizes the 
computed antenna performance at 0.22 f . The antenna efficiency for 
the TE U mode is relatively high particularly in view of the large 
phase variations of the subreflector illumination. The far field radia- 
tion patterns for both the TE n and TM i modes show good charac- 
teristics. The primary disadvantages, however, are the high noise 
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Fig. 14 — Horn radiation pattern at distance S (TEn mode, freq. = /„). 
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Fig. 15 — Far field radiation pattern (TE n mode, freq. = /<>). 

temperatures for both modes, owing to the power loss at the subre- 
flector. 

VII. SUMMARY AND CONCLUSIONS 

Computer programs have been developed for computing the radia- 
tion characteristics of Cassegrain antennas and for plotting of the 
computed radiation patterns. The method is applicable to symmetrical 
Cassegrain antennas and provides the means of their design for nearly 
optimum performance. 

A Cassegrain antenna with a near field conical feed horn has been 
investigated for different mode excitations and over a wide frequency 
range. For large antennas (over 200 wavelength main reflector diam- 
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Fig. 16 — Horn radiation pattern at distance S (TEu mode, freq. = 0.22 x /»). 
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Fig. 17 — Far field radiation pattern (TEu mode, freq. = 0.22 x /.). 

eter) this type of feed can be used over a 50 per cent bandwidth 
with only small variations in the over-all antenna characteristics, 
except for the predictable increase in the antenna gain with frequency. 

The computed antenna characteristics for the combined TEn and 
TMn mode excitations show that the advantages of the combined 
excitation are: (i) lower far field E-plane sidelobes, (n) 0.2 dB 
higher antenna gain, and (Hi) 10°K lower antenna noise temperature. 

At the design frequency for TEn mode excitation, the computed 
antenna efficiency is 70 per cent and the noise temperature near hori- 
zon 15.5°K. With a design modification it should be possible to 
reduce the noise temperature by a few degrees without affecting the 
antenna gain. 

For the TM i mode, the calculated antenna gain and noise tem- 
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Fig. 18 — Horn radiation pattern at distance S (TMm mode, freq. = 0.22 X /»)■ 
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Fig. 19 — Far field radiation pattern (TEn mode, freq. = 0.22 X /«)• 

perature at the design frequency are superior to those obtainable by 
using a conventional feed. The sidetobe levels of the far field radia- 
tion pattern are perhaps a few dB higher. 

At a frequency 4.5 times below the design frequency the calculated 
antenna efficiencies for the TE U and TM i modes are relatively high. 
However, the power loss at the subreflector gives rise to appreciable 
noise temperatures near the horizon. 

APPENDIX A 

Formulations Used in Computing the Characteristics of Cassegrain 
Antennas 



A.i Horn Radiation Patterns 

The horn radiation patterns have been computed by using the Kirch- 
hoff approximation to the aperture radiation field. With this approxi- 
mation the electric field E p at distances of at least a few wavelengths 
from the aperture is \ 

Ep = fr // [Ea(1 + ln ' u " E -' lft(ln + 1r)] St^ s (10) 



where 



S = horn aperture area 

k X 

X = wavelength 
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and 1„ and l w are unit vectors in the normal, and in the R direction as 
shown in Fig. 20. E u is the electric field in the horn aperture, assumed to 
be the same as for circular waveguide modes but with spherical wave 
fronts 

Since it lias been shown by actual computations 4 that the primary 
contributions to the horn radiation patterns are due to the first terms 
of (10). in the computations the term E-1 B (1„ + 1 R ) has been neglected. 

Because of the periodicity of the e~' k R/R with respect to the azimuth 
coordinates <?' and <f , , it is sufficient to evaluate (10) in a discrete 
number of ip t planes, the number I eing equal to number of Fourier 
components of the aperture field in c . In particular it has been shown 1 
that for TE n and TM U mode excitations it is sufficient to evaluate one 
rectilinear z or y component of (10) in the two principal planes <p v — 
and <p x = 7r/2. Similarly for TM , mode excitation only one rectilinear 
component of E„ in one plane need* to be evaluated. 

The integrals which are evaluated I'm the TE xl , TM X \ , and TM 01 
modes are: 

E„ = k £ I" f'z-^Q- + i"- 1 *)*' 11 o'de'dv' (ii) 

with the aperture fields for the different modes given by: 
TE n mode 



CM* 



= J.\k 



- J 2 [ k TEl , — ) cos 2tp' 



(12) 



P(r,', ff,',*,' 
R P(r„*,,0i) 




Fig. 20 — Coordinates for horn radiation pattern computation. 



916 THE BELL SYSTEM TECHNICAL JOURNAL, JULY-AUGUST 1068 

with 

J{(itB„) = 0. (13) 

TM,, mode 

(# fl y)TM„ - J.\kTM„ 5/ + J2 ( A "™" 7J C ° S V (14) 



with 



TM 01 mode 



/iGfew.,) = o. (15) 



with 

/.(*,«„) = 0. (17) 

J„ = Bessel functions of order n. 

a = horn flare angle. 

The integration with respect to / has been eliminated by approxi- 
mating the integrals, I nh , given by: 

-2T -jkR 

I nh = J ~ (1 + ln-U) cosnfo - <p') d(jpx ~ <P')> (18) 

The approximations used are modifications of the previously derived 
approximations, I\ , to the integrals (18) with l n - 1 R = 0. The modifica- 
tions consist of including the values of l n -l« at the stationary phase 
points since it has been shown that the previously derived approxima- 
tions, I\ , reduce to those obtained by the method of stationary phase, 
and that I\ can be separated into terms which correspond to the sta- 
tionary phase terms. It is subsequently shown that R can be expressed 
as: 



R = r Jl-y cos ( ft - <p') (19) 

where r and u are functions which are independent of ip and <p', on this 
basis, the first order approximations to (18), V nh are: 
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-jkdlt-u) 



/-'/, = 7rf|L R (1 + ln-l«„) + — ^ (1 + 1,,-W y n (ku) 

- j[ Rg (i + i«-i«J - - Ri (i + 1»-UJJ^W| (20) 

with 

R ° = 'i 1 ~ vf (21) 

Rt = /(l + yf (22) 

and l* o and l Ri are the unit vectors at (<px — <p') equal to zero and w 
respectively. 

By using the approximation (20), one integration is eliminated in (11 ) 
and the radiation patterns for the different modes are computed from 
the following integrals : 
TE U Mode 

(#Jte„ = *~ [ [./„(a- T k„ 5)/L ± ./ 2 (*™„ 5)^] «'» *' W 

(23) 

where the minus signs give the radiation pattern in the plane ^ = 
and the plus sign the radiation pattern in the plane f t = ir/2. 

TMn mode 

The integral is analogous to (23 1 . 

TMot mode 

(EJ TMoi = ^ /" ./,(a- TMo1 £)l\ k sin d' dd'. (24) 

Referring to Fig. 20, 

R = (I 2 + r 2 + p 2 + 2rjp cos 0, - 2lp cos 0' - 2r,Z cos tO* (25) 

with 

cos 7j = sin X sin 0' cos (p' — ^,) + cos 0! cos 0' (26) 

and 

i i _ 2 cos g/ ~ l + r ' cos 7 1 / Q7 x 

lnl« - fl (27) 
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Hence a comparison of (25) with (19) shows that 
r = (I 2 + r? + p 2 + 2r,p cos 0, - 2lp cos 6' - 2r x l cos 0, cos 0')* (28) 
and 



u = 



rjs'm 0, sin 6' 



(29) 



The integrals for the difference modes have been computed at two 
values of r x \ (i) at the subreflector surface to obtain the subreflector 
illumination, and (ii) at a constant distance corresponding to short- 
est distance, s, from the subreflector to the horn aperture. The latter 
was performed to obtain the horn radiation pattern in a form which 
is readily measurable and convenient for subsequent computation of 
the horn power radiation pattern used for determining the power loss 
at the subreflector. 

For the first computation referring to Fig. 21 



r, = 



c (1 - /3 2 ) 
2 cos 0,-/3 



(30) 



with 



'-.*■ 




Fig. 21 — Subreflector coordinates. 
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For determining the subreflector illumination it is preferable to 
obtain the illumination in terms of subreflector coordinate 2 . The 
relationship between the coordinates 6 X and 8 2 is: 

(1 + /3 2 ) cos d 2 + 2/3 
C0S 6l 1+^ + 2/3 cos 2 (31) 

The second computation 



c 



r, = \ (1 + 0) (32) 

and the integration is performed as a function of 6 X . 

A comparison has been made between some radiation patterns 
computed by single and double integration. Good agreement was 
obtained. 

The electric field in the spherical 6 X and <p x coordinates can be ex- 
pressed in terms of the radiation patterns in the principal planes 
ipx = 0, it/2. 

ForTE„andTM n modes 



E„ = U t E„\Z) sin Vl + l Pl E vy (0) cos <p, (33) 

The TMoi mode has only a X component given by (24) . 

a.2 Horn Power Radiation Patterns 

The horn power radiation pattern, P, r , is computed from the fol- 
lowing integral: 

a -2» -9, 

P " = 2~ J J n EpE * S ' m 6 dd d,f - < 34 ) 

■q = free space intrinsic impedance. 

The total power is obtained by extending the range of t to the region 
where E p -E* has a negligible value. 

The total radiated power can also be obtained from the assumed 
fields at the aperture (12) through (17). On this basis the total 
power for the different modes is approximately 

TE n mode 



TM n mode 



*V L A' TEll J 



(35) 



r, 2ir(la) r o,. 

P ", = ~ Jp^TM,,)' (36) 
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TMoi mode 

Plli ^^fj^ kTM J. (37) 

The computations performed for the different modes and fre- 
quencies by using (34) are in agreement with (35) through (37) 
within 1.8 per cent, with the power computed by using (34) giving 
larger values for all modes. This is at least partially caused by the 
higher values that the approximations I,, 1 give compared with those 
obtained by precise numerical integration. 2 

A.3 Subreflector Radiation Patterns 

The subreflector radiation patterns have been computed by using 
the surface integral relating the radiated fields and the current dis- 
tribution over a surface. 7 For the distance of at least a few wave- 
lengths from the reflector the radiated electric field with reference 
to Fig. 21 is: 

E - = x // lw - X (J X 1sb) ' IT d8a (38) 

where 

J = surface current density 
S, = subreflector area. 

To evaluate (38) it has been assumed that the reflector is locally 
plane. With this assumption the current density is directly related to 
the incident electric field. To simplify the computations 1*. was re- 
placed by l r3 . (A test computation of the subreflector radiation pattern 
of the open Cassegrain antenna showed that using l r , instead of 1 R , 
results in a negligible difference.) 

For a hyperboloid reflector, the relations between the incident 
electric field and the current density have been derived. 1 By using 
the approximation (20) the integration with respect to <p 2 can be 
eliminated. On this basis the radiated fields of the subreflector result- 
ing from the incident fields of a TE n mode (23), are given by: 

In the plane y> 2 = 0, 
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In the plane f> = n/2, 

7T 

2> 



EAl 



_L 
2X 



I e \L E * M + MU + cos £ > "" * 

Tet aa /■ M (1 + /3 cos 0,)~L 
- [KM - B„(-J - + P cQs ^ j/; cos » 

+ 2, --(i) ^b; /: sin ')* sin e - de > • (40) 

For the TMn mode, the fields are analogous to (39) and (40) 
with the corresponding TM ai illuminated functions E p (0) and E p (tt/2). 
For the TM i mode, 

[£.(0)] TMol = -^ fj E P M e- ik " 

i „ (1 4- cos 2 ) ,, . /3sin 0sin 0, rl ~| ■> • „ ,„ #-,% 

• cos 5 o _i_ — a 7 + Vt — -r 7 » H 8m ** d$ * ( 41 ) 

L /3 + cos 2 /3 + cos 2 J 

where E p (0) is given by (24) . 

In the above r 2 is the equation of the subreflector in the coordinates 
shown in Fig. 21, and is given by 

c (1 - /3 a ) m 
12 2/3 4- cos 02 (4J} 

I\ are the first order approximations to the integration with respect 
to <p 2 . They are given by (20) with 1„- \ Ro and l„-l Rl set equal to zero. 
The other parameters are: 

Ri - (/•■"' + ri - 2rr 2 cos 72)*. (43) 

cos 7 2 = sin sin 2 cos (<p — <p 2 ) + cos cos 2 . (44) 

/To sin sin 0o . 

U2 ~ (r 2 + r 2 2 - 2/r 2 cos cos 2 ) J ' (45) 

The subreflector radiation patterns have been computed at two 
distances: at the paraboloid surface 

r - rrk-e <« 

where / = focal length of the main reflector, to obtain the main 
reflector illumination, and at 

r = / (47) 

to determine the power loss at the main reflector. 



922 THE BELL SYSTEM TECHNICAL JOURNAL, JULY-AUGUST 1908 

The if dependence of the subreflector radiation patterns are ob- 
tained in terms of radiation patterns in the principal planes and are: 

For TEu or TM„ modes 

[E.] T k,,.tm„ = 1^.(0) cos<p + !•£-(§) sin *»■ ( 48 ) 

For the TM , mode 

[E.] TM „, = 1.^.(0). (49) 

The power radiation pattern of the subreflector is computed by using 
the integral (34). 

a.4 Aperture Gain and Efficiency 

The aperture gain and efficiency are computed by projecting the 
incident field on the main reflector aperture. The fields in the recti- 
linear x, y components are related to 6, ip components by the ex- 
pressions 

E T , = —1 X (E.„ cos <p — E t9 cos <p) 

- l„(E,i sin <p 2 + E. r cos <p 2 ). (50) 
For TEn and TM 1X mode excitations polarized in the y direction, 
the gain on axis, Om , is 



G u - x , 



["" \ E, \M + E '®J sin 6 d ° 



C \\ e $ 2+ \ e " (o) \j s ' mdde 



(51) 



where r is the equation for the paraboloid (46) . 

The aperture efficiency, g, is obtained from the relation 



9 = 



Q 

G 



(52) 



with 

G n = 



f" 4tt/ sin 6 m T , r)3 , 

LX(1 + cos O J 



The maximum antenna gain for the TM i mode is determined by 
normalizing the amplitude of the TM i mode electric field at the 
horn aperture with respect to the amplitudes of the TE n or TEn 
and TMn modes for the same power input, by using (35 through 37). 
The gain for the TM i mode is then related to the gain for the TE n 
mode on axis, referred to the maximum of its pattern. 
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A. 5 Far Field Radiation Patterns 

The far field patterns are computed from the projected field on 
the aperture, using the relation 



• / / E r , exp [jkr sin sin 0„ cos (<p — <p a )]r 2 sin dd dtp (54) 

where r„, 0„, and <p a are the coordinates of the far field observation 
point. 

Because of the antenna symmetry the integration with respect to 
if can be the readily performed. The resulting integration with respect 
to is for TE„ or TM U modes 

± [j?,(|) - £\(0) J./ 2 (^'sin sin a )}r 2 sin f/0 (55) 

where =fc signs correspond to the patterns in the H(ip a = 0) or E(<p a = ir/2) 
planes. 
Similarly for the TM , mode 

(E/) TM o, = U. y e -^r~ / A'»(0)./,(^'sin sin o )r 2 sin d0 (56) 

where Z? 8 is the main reflector illumination (48) and (49) for the 
different modes in the planes <p = and f = tt/2 . 

APPENDIX B 

Program for Computing the Characteristics of Cassegrain Antennas 
and for Graphic Display of Radiation Fields 

The package consists of two programs: a program to compute the 
horn radiation patterns, the subreflector radiation patterns, the far 
field radiation patterns, and other characteristics of Cassegrain an- 
tennas described in Appendix A, and a program to scale, label, and 
plot the radiation fields. The two programs are linked by intermediate 
storage of computed results and control variables on tape, at the 
conclusion of Part 1 execution. 

B.i Computation Program 

Figure 22 is a logic diagram of the program. The following conven- 
tion is used throughout the logic diagram. Square-bracketed symbols 
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Istart7 



READ NAMELIST 
PDATA 




READ NAMELIST 

HDATA 

I 

(DO 50 KK = i,2j 



THS[« 2 ]=o.o 
THC [0|]=o.o 



CALL SIM, TO CALC. 
EHS(N)[(E py )]0=o 

EES(N)[(E py )]0=w/2 

ETMS (N) [(E py )] TM 01 
AT SUB RE ELECTOR 



CALL SIM, TO CALC. 

EHC(N)[(E p y)]«=o 

EEC(N)[(E py )]^=^/2 

ETMC(N) [(E py )]TM 01 
AT DISTANCE S 
I 



( DO 50 N = i,L_) 

1 THH[g]=0.0 1 

(DO 60 J = i , nt) 
I 



CALL JOI2.CALC. 
APERTURE FIELDS 

( E ay) T E M '( E ay)TM (l 

( E ay) T M i 



CALC. 

RL[r,tt g )] 

CTHi[cos{0,(0 2 )}] 
I 



CALC. 

RL[r,(fl 2 )]= const 
CTHi[cosW0] 



CALC 

RN[r], X[u] 
R0[R ], Ri[R,] 
APL[]n,lR ] 

AMI [t n ,l Rl ] 

I ~ 



CALL 1012, CALC. 

IO[I'oh],Il[l',h] 

I2[l 2h ] 




EVALUATE INTEGRAND 
AT SUBREFLECTOR 

Si (J) H-planeI TE„ 

S2(J)E-planeJtm„ 

S3(J)TMoiMODE 



evaluate integrand 

at distance s 

S4(J)H-plane]TE„ 

S5(J) E-PLANEjTM,, 

S6(J) TM ot MODE 



THH=THH+A 



THC = THC+A 3 | 




I READ NAMELIST SUBDAT| 
I 



CALL QUAD, TO CALC. 

EHSU (NTS)=QUAD[EHS CD] 

EESU (NTS) = QUAD[EES (L)] 

ETMSU (NTS)=QUAD[ETMS (U 



(DO 500 MM =1,2) 

, I 

TSP[9]=o.o | 
I 
(DO 500 N = i,LS ) 



TSU[fl 2 ]=0.0 



(DO 510 K = 1,NTS) 
... -^ N 
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£ 



calc. RHC[r 2 ] 

X5P[U 2 ] ETC. 



CALL 101 2. CALC. 10 S [1'^ 
IlS[l',]l2S[l' 2 ] WITH 
ln-tRo =, n-'Ri = 




EVALUATE INTEGRAND 
AT PARABLOID SURFACE 



SHP(K)H- PLANE] 
SEP (K)E- PLANE 



TE„ 
TM« 



STMP(K)TM 01 MODE 



hVALUATE INTEGRAND 
AT DISTANCE f 

SHPCCWH-planeVtE,, 

SEPC(K)E-planeJTM„ 
STMPC(K)TM , MODE 




CALL SIM, CALC. 
EHPCN)[(E s )]^ = o 

EEP(N)[(E s )]0=ir/2 

etmpcn)[Ce s )]tm i 

AT PARABLOID SURFACE 



CALL SIM, CALC. 
EHPC(N)[(E S )]0=O 

EEPC(N)[(E 5 )] t>=ir/2 

etmpc(n)[(e s )]tm 01 

AT DISTANCE f 



TSP = TSP-(-A* 



TSP = TSP+A 6 




1 



READ NAMELIST 
FFDAT 

I 



CALL QUAD, TO CALC. 

EHPF(NTP) = QUAD[EHP(LS)] 

EEPFCNTP)=QUAD[EEP(LS)] 

ETMPF(NTP) =QUAD[ETMP(LS)] 



TPF[9 2 ]=0.0 



( DO 700 M= I, LP ) 

TPD[g] = 0.0 
( DO 7io j==i,NTP) 



EVALUATE INTEGRAND 

shpf(j)h-plane1te„ 

s epfcj)e- plane tm,| 
sthpf(j) tm . mode 



(7io)- JTPD = TPD + A 7 



CALL SIM, CALC. 

EHFFCM)[(EfgJ*a=o 

EEFF(M)[(E fy )]0 a =ir/2 
ETMFF(M)[(E f )]TM 0l 




Fig. 22 — Logic diagram for field calculation program. 
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follow the notation of Appendix A, while the preceding alphanumeric 
symbols are the Fortran IV source program names. 

Four groups of input data are required, designated by the namelist 
names pdata, hdata, subdat, and ffdat. Although some of the fol- 
lowing data is redundant, the formats are designed for convenience 
and precaution. 

PDATA 

ns — Number of sets of field patterns to be computed, that is, 

NS = 2 Horn radiation patterns only 
NS = 4 Horn and subreflector patterns 
NS = 5 The above plus far field patterns 

nps — Number of patterns per set, that is, 

NPS = 2 E & H planes only 

NPS = 3 The above plus TM 01 mode 

ite — Control bit for plotting program (see Section 2) 

HDATA 

hl — Horn length 

hlamd — Horn length normalized with respect to design wavelength 

fq — Frequency normalized with respect to design frequency 

c — Distance between foci of hyperboloid 

beta — Defined by equation (30) 

pl — Location of phase center normalized with respect to horn 

length 
alpha — Horn flare angle 
l — Number of points at which horn radiation pattern will be 

evaluated 
nt — Number of points at which integrand will be calculated for 

evaluation of integrals in equations (23) and (24). 
deg — Angular increment (in degrees) for obtaining subreflector 

illumination 
ang — Angular increment (in degrees) for obtaining horn radiation 

patterns 
tm„ — Complex constant which determines TM U mode to TE„ 

mode ratio 
up — Control bit for plotting program (See Section B.2.) 
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SUBDAT 

f — Focal length of paraboloid 

Note: The dimensional unit for f, c, hl must be the same. 
flamd — Focal length normalized with respect to design wavelength 
clamd — C normalized with respect to design wavelength 
ls — Number of points at which subreflector radiation pattern 

will be evaluated 
nts — Number of points at which integrand will be calculated for 

evaluation of integrals in equations (39), (40), and (41) 
gama — Main reflector illumination angle (degrees) 
degp — Angular increment (in degrees) for obtaining main reflector 

illumination 
angp — Angular increment (in degrees) for obtaining subreflector 

radiation patterns 
incs — (AT + 1) where N is the number of points to be interpolated 

between previously computed subreflector illumination points 
lb — Number of points which are not included in integral (51) 

because of subreflector blocking 
i2P — Control bit for plotting program (See Section B.2.) 

FFDAT 

ntp — Number of points at which integrand will be calculated for 

evaluation of integrals in equations (55) and (56) 
gamab — Angular portion (in degrees) of main reflector blocked by 

subreflector 
ntpb — Number of points which are not included in integrals in 

equations (55) or (56) owing to subreflector blocking after 

interpolation of main reflector illumination 
lp — Number of points at which far field radiation patterns will 

be evaluated 
incp — (N + 1) where N is the number of points to be interpolated 

between previously computed main reflector illumination 

points 
degff — Angular increment (in degrees) at which far field radiation 

pattern will be evaluated 
I3P — Control bit for plotting program (See Section B.2.) 

The following subprograms must be included in the deck before 
execution. 
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prad — Computes the power radiation patterns in accordance with 
equation (34) 

jo 12 — A Bessel function subroutine developed by J. Alan Cochran 
and P. A. Alsberg. The present version includes two sub- 
sidiary subroutines : 8 

dphase — Uses phase-amplitude method for large values 

of argument 
jlow — Uses downward recursion technique for small 

values of argument 

1012 — Calculates the first order approximations V„ h (20) to the 

integrals I„ h (18) 
quad — A quadratic interpolation scheme for complex arrays 
interp — A subroutine called by quad 
proc — A subroutine to format and print output data 
sim — A complex Simpson's rule integration routine. Will accept 

an even or odd length array with negligible variation in 

accuracy 
gain — Computes the antenna gain and aperture efficiency as defined 

by equations (51), (52), and (53) 
tr — A special purpose of Simpson's rule integration to evaluate the 

integrals in equation (51). This function subprogram is called 

only by the gain subroutine 

The program requires approximately (52,660) 8 or (22,000) i words 
of storage. A representative execution time for both modes (TEn 
and TMn combined, and TM i) at the design frequency is 14 min- 
utes; the same calculations at 0.22 times the design frequency, where 
a smaller number of integration points is required, takes approxi- 
mately 5 minutes. 

B.2 Plotting Program 

A logic diagram of the program is presented in Fig. 23. 

All input data required by the plotting program has been stored 
on tape by the previous program. 

The plotting control bits, referred to in Section 1, have the fol- 
lowing meaning: if field calculations are to be made in a combined 
TE n and TMn mode— that is, input data TM n ¥= (0.0, 0.0)— the 
control bit ite in namelist pdata must be set equal to 2. If the field 
calculations are to be made in the TE n mode alone — that is, TM U = 
(0.0, 0.0)— the control bit should be set equal to 1. Therefore calcu- 
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lations for the three sets of radiation patterns, horn, subreflector, 
and far field, will generally be made in two modes, a combined TEu 
and TMn mode, and the TM fll mode, where it is understood that the 
combined mode may be the pure TE n mode, if TM n = (0.0, 0.0) and 
ite = 1. 

For the combined mode the radiation fields will be evaluated in 
both E and H planes. E- and H-plane data are plotted together for 
ease of comparison. However, in some cases (particularly for certain 
far field patterns) a rapidly varying phase plot superimposed on a 
rapidly varying amplitude plot may result in an unclear graph. For 
this reason the control bits up, i2p, and I3P are introduced, up con- 
trols horn radiation pattern plotting, I2P controls subreflector plotting, 
and I3P, far field plotting. If the control bit for a particular field is 
set equal to 0, two plots will be generated, that is, 

Combined Mode 

(i) E-plane and H-plane amplitude and phase 

TMoi Mode 

(u) Phase and amplitude 

However, if the control bit is set equal to 1, four plots will be 
generated : 

Combined Mode 

(i) E-plane amplitude and H-plane amplitude 
(ii) E-plane phase and H-plane phase 

TMoi Mode 

(Hi) Amplitude 
(iv) Phase 

Vertical scales are restricted to allow only 20 divisions, therefore a 
preferred set of increments for the various scales has been selected. 
The allowed increments in dB for the amplitude scale, stored in 
array adb(i), are: 

0.5, 1.0, 1.5, 2.0, 2.5, 3.0, 4.0; 
for the phase scale, in degrees in array aph (i) : 

1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 8.0, 10.0, 12.0, 15.0, 18.0; 
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vrt/ 



(DO 73 J = i, NS) 



(DO 72 K = 1,NPS) 



K=i: ER(I) = ERi (I); EI(I) = EI(I) 
K = 2: ER(I) = ER2(I); EI(I)=El2(I) 
K=3: ER(I)=ER3(I); EI(I) = El3(I) 



amp (1 ) = [[ER (i)] 2 + [EI (I)] 2 ] ,/z 




FIND N 

amp(N) = max[amp(I)] 

1=1, K2 



EN = AMP(N) 
EKN) 



PN = TAN 



ER(N) 



EN = AMP0) 
EIO) 



PN = TAN 



ER(D 



INTERPOLATE 
REAL PART ER(I) 



INTERPOLATE 
IMAG PART EI (I) 



CALC. RELATIVE AMP 
dB (I)= 20.0 X LOG ^gp 



CALC. RELATIVE PHASE 

RP(I)=TAN-' l^^l-PN 



Ler(i)J 



CONSTRUCT AMP SCALE 

YU2 max {dB(D} 
YSi^ min {dB (I)} 



CONSTRUCT ANGLE OFF 
AXIS SCALE 




SYLi=YLi 
SYSI=YS1 
SYL2=YL2 
SYS2 = YS2 
SDB(I)=DB(I) 
SRP(I)=RP(D 




N ^PT> Y 



COMPLETE TITLE 
FOR PLOT 



COMPLETE TITLE 
FOR AMP PLOT 



CALL FRAME CALL FRAME 



CALL PLOT 2 



PLOT DB(I) 



PLOT 

DB(I), RPU) 



COMPLETE TITLE 
FOR PHASE PLOT 



CALL FRAME 



PLOT RP(I) 



CONSTRUCT PHASE SCALE 

Yl_2£ MAX {RP(D} 

YS2=-YI_2 
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Fig. 23 — Logic diagram for field plotting program. 
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for the angle-off-axis scale, in degrees, stored in array as(i): 

6.0, 5.0, 4.0, 3.0, 2.5, 2.0, 1.5, 1.0, 0.75, 0.5, 0.4, 0.3, 0.2. 
The following subroutines must be included in the deck before 
execution : 

plot 2 — A subroutine to generate a grid with two independently 

labeled ordinates sharing a common abscissa 
minmax — A subroutine to select the algebraicly largest or smallest 

entry in an array and specify its index 
interp — A quadratic interpolation scheme for real arrays 
filter — Adjusts plotting data for phase variations in the vicinity 

of ±180° 
label 2 — A modified version of the microfilm subroutine label. 

Called only by plot 2. 
The program requires approximately (51,536) 8 or (22,000) w words 
of storage, and about 0.5 minute execution time for all twelve plots. 
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